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Photons emitted from the electromagnetic fields of relativistic heavy ions 
' can fluctuate into quark anti-quark pairs and scatter from a target nucleus, 

m 

' emerging as vector mesons. These coherent interactions are identifiable by 

a^ . 

' final states consisting of the two nuclei and a vector meson with a small 

transverse momentum. The emitters and targets can switch roles, and the two 
possibilities are indistinguishable, so interference may occur. Vector mesons 
^ • are negative parity so the amplitudes have opposite signs. When the meson 

transverse wavelength is larger than the impact parameter, the interference 
j_j ■ is large and destructive. 



p. 



The short-lived vector mesons decay before amplitudes from the two 
sources can overlap, and so cannot interfere directly. However, the decay 
products are emitted in an entangled state, and the interference depends on 
observing the complete final state. The non-local wave function is an example 
of the Einstein-Podolsky-Rosen paradox. 
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In relativistic heavy ion collisions, vector mesons are copiously produced via photon- 
Pomeron fusion [|T|. The photon and Pomeron both couple coherently to their emitting 
nuclei, giving these reactions a distinctive signature, with the final state consisting of the 
two nuclei, a vector meson with low perpendicular momentum {p±) and nothing else. The 
latter requirement restricts these interactions to peripheral collisions, usually with impact 
parameter b > 2Ra, Ra being the nuclear radius. At heavy ion colliders, vector mesons with 
masses up to about 2'yhc/RA, where 7 is the Lorentz boost of each beam are produced. 

The electromagnetic field has a long range, while the nuclear (Pomeron) field has a short 
range compared to the size of the nucleus. So, vector meson production occurs in the region 
occupied by the Pomeron emitting ('target') nucleus. Since the production is coherent over 
the entire target nucleus, it is a fairly good approximation to treat the meson production 
regions as two point sources, one at the center of each nucleus. The situation is similar 
to that in a two-source interferometer, albeit with unstable particles. A parity inversion 
switches the emitter and target. Because of the vector meson negative parity, the two 
emission amplitudes have opposite sign. For perpendicular momentum p±_ <^ 1/6, b being 
the impact parameter, the interference is destructive. We calculate here the magnitude of 
the interference and discuss the implications of the short vector meson lifetime, which causes 
the mesons to decay before wave functions from the two sources can overlap. 

The cross sections for meson production can be calculated by convoluting the Weizsacker- 
Williams virtual photon spectrum with the photonuclear interaction cross section. The 
photonuclear cross section is determined from data on 7p interactions 0, using a Glauber 
formalism to find the cross section for — > VA, where A is the target nucleus, and V a 
vector meson. In this approach, the photons are considered to fluctuate to virtual vector 
mesons, which then elastically scatter from the target nucleus, emerging as real mesons. 

Earlier calculations of production cross sections summed the cross sections from the two 
nuclei and found very large cross sections [0. Exclusive p° production was about 10% of 
the hadronic cross section for gold on gold collisions at a per nucleon center of mass energy 
\^Snn = 200 GeV as will be studied at the Relativistic Heavy Ion Collider (RHIC), rising 
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to 30% for the a/Stvtv = 5.5 TeV lead on lead collisions at the Large Hadron Collider (LHC). 
The uj and (f) production are about an order of magnitude smaller, with J/ip cross sections 
of 0.3 mb (gold at RHIC) up to 32 mb (Pb at LHC). These cross sections correspond to 
production of hundreds of p° per second at RHIC, and hundreds of thousand of p° per second 
with calcium beams at the LHC. The J/ip rates range from 0.06 Hz with gold at RHIC up 
to 780 Hz with calcium at LHC. These rates are comparable to those found at current and 
future meson factories. 

Because of the small 6, for low final states it is impossible to tell which nucleus emitted 
the photon, and which elastically scattered the meson. So, the amplitudes are combined. 
Because vector mesons are negative parity, the two amplitudes subtract, rather than add, 
leading to destructive interference. In different terms, for p±h < 1, a system of two identical 
nuclei has zero dipole moment, so radiation of vector particles is impossible. This situation 
is analogous to bremsstrahlung by identical particles |^ , where there can be almost complete 
destructive interference. However, in contrast to bremsstrahlung photons, vector mesons are 
short lived, and they don't live long enough to travel the distance h. So, the vector mesons 
generally decay before their wave functions can overlap. However, their decay products can 
interfere. 

The amplitude for vector meson production is A{y^ p±, 6)e*'^*^^\ with A(y, p±, b) the magni- 
tude and (f){y) the phase at rapidity y. We assume that A{y,p±, b)e^'^^y^ is symmetric to rota- 
tions around the beam direction. The photon energy and perpendicular momentum are k and 
k±, and the final meson momentum p. All values are in the center of mass frame, which cor- 
responds to the laboratory frame for a heavy ion collider. To eliminate the directional ambi- 
guity, we adopt the convention that, for y > 0, the photon energy is higher than the Pomeron 
energy. Then, A{y > 0) < A{y < 0) and y = ln(2A;/My). These amplitudes are determined 
from the cross sections calculated in Ref. using A{y,p±,b) = o-A+A-^A+A+vijJi P±,b), 
for a single production direction. 

The final state perpendicular momentum is the sum of the photon and Pomeron per- 
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pendicular momentum; the spectrum is a convolution of the two sources. The photon 
perpendicular momentum can be found with the equivalent photon approximation 



(Pk^dk 7r^{kl + ky-f^y 

The nuclear form factor, 

1 



(1) 



^^^^ = -J^ 



sin(gi?^) — qRACOs{qRA) 



1 + a2g2 



(2) 



the convolution of a hard sphere and a Yukawa potential with range a = 0.7 fm , is an 
excellent fit to a Woods-Saxon density distribution Here, po is the nuclear density and 
A the atomic number. For a given k, k± is independent of b, with d^N^{k, k±)/d'^k±dk rising 
from at /cx = to a maximum at k± = k/'y and then dropping as k± rises further. 

Since the Pomeron range is short, its perpendicular momentum spectrum is determined 
by Eq. (0). Figure 1 shows the photon and Pomeron perpendicular momentum spectra, 
along with their convolution, for y = 0, and y = —2. Diffractive fringes are visible in the 
Pomeron spectrum; since the photon contribution to p± is usually small, these fringes also 
appear in the final p± spectrum. 

The production phase depends on the process. Since the soft Pomeron represents the 
absorptive part of the cross section, the amplitude for photon-Pomeron fusion should be 
largely imaginary, with a small real part because the cross section rises slowly with photon 
energy Parameterizing aljN —>■ VN) ~ W^, where W is the 7-nucleon center of mass 
energy, (f){y) = tan^-'^(7re/4), independent of y 0. For the light mesons, e = 0.22 so ~ 10°. 
For the J/ip, cr{'~fN VN) ~ W'^'^, inconsistent with the soft Pomeron; the steep rise may 
be due to threshold behavior or signal the breakdown of the soft Pomeron model |^; in 
either case, the J/ ip phase angle must be treated with caution. 

For the and cj, photon-meson (primarily the /o 0) fusion also contributes to the cross 
section, which may be parameterized as XW^ + YW~^, where 77 > 1, so the photon-meson 
contribution decreases as W rises. Besides increasing the cross section, this causes (j){y) to 
vary with y. Because the phase has only been measured as an average over photon energy 



P], we will treat (j){y) as a constant. Of course, at ?/ = 0, the energies are equal; as \y\ rises, 
the phase difference could become significant. 

For two nuclei at points xl and X2, the amplitude Aq for observing a vector meson at a 
distant point xq is found by approximating the vector mesons by plane waves. The meson 
momentum p is determined by p± and k via p|| = /c — rriy/Ak. 

Ao{xo, p, h) = A{p^, y, 6)e*[</'(3/)+P-(^l-^<')l _ A{p^, -y, 6)e^['^(-?')+P-(^^2-xo)] (3) 

with the negative sign because of the negative parity. Because of their short lifetimes, the 
vector mesons will decay before reaching xq- In fact, they will generally decay before their 
wave functions can overlap. However, their decay products can interfere; with p now being 
the sum of the product momenta. The resulting virtual interference pattern depends only on 
p- b. The decay is incorporated by multiplying A{p^, y, h) by the decay amplitude, including 
the branching ratios, lifetime and angular dependence. Since the angular distributions are 
the same for the two production directions, the decay amplitudes are the same for both 
directions, and does not affect the observed cross sections. 

Since the decay products fly off in different directions at relativistic velocities, the pro- 
duction amplitudes from the two sources can only interfere if the final state wave function 
is nonlocal P]. For example, in J/ip ^ e+e^, the electron and positron are nearly back-to- 
back, and the e"*" (or e~) amplitudes from the two sources cannot overlap each other until 
they are a good distance from their parent J /ip. By the time the (and e~) from the two 
sources overlap, the and e~ are well separated, and any interference involving both the 

and e~ requires a non-local wave function. 

The final state wave function is entangled, with the form exp(z/ci-xi) exp(iA;2"3^i)— exp(iA;2- 
X2) exp{ik2 ■ X2) where ki and k2 are the individual decay product momenta (for a 2 body 
decay), and cannot be factored into individual particle wave functions. This interference 
is visible by observing the complete final state, but not by examining the individual decay 
products. This non-locality is an example of the Einstein- Podolsky- Rosen paradox | |10[| . 
However, the correlation depends on p±, which is a continuous variable, rather than a 
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discrete variable like spin, as is commonly considered. 

The individual decay products have much higher momenta than their parent meson, and 
in principle could be used to reconstruct the meson decay point, and hence determine which 
nucleus emitted the meson. However, the position measurement would obscure their momen- 
tum enough to eliminate the interference pattern, as occurs with a two-slit interferometer 

n . 

In the center of mass frame, production is nearly simultaneous (within At = b/cy), so 
time drops out. Since A is defined to be real, 

^{P±, y, b) = y, b) + -y, b) - 2A(p^, y, b)A{px, -y, b) cos[0(?/) - (j){-y) +p-b] 

(4) 

At midrapidity, y = 0, the two source contributions are equal and the observed cross section 
is 

a{p^,y = 0,6) = 2A\p^,y = 0,6)(1 - cos[p- b]). (5) 

For a given 6, a oscillates with period Ap_|_ = h/b. When p ■ b <^ 1, the interference is 
destructive, and there is little emission. 

Since b is in principle measurable by examining the outgoing ions, it is an observable, and 
we integrate the cross section (not amplitude) over all possible b to get the total production. 
Figure 2 shows (PN/dp\ at y = for gold and silicon beams at RHIC and calcium at LHC, 
with and without interference. 

In theory, one could choose to measure the perpendicular momentum transfer from the 
two nuclei instead of b. Then, for roughly p± > 10 MeV/c, the larger momentum transfer 
will usually come from the target. This determination will break the symmetry, reducing 
the magnitude of the interference. Of course, near pj_ = 0, the photon and Pomeron per- 
pendicular momenta have very similar magnitudes, so the interference will remain. Because 
the initial nuclear momenta are not well known, this possibility seems at best theoretical. 

With gold beams at RHIC, vector meson production occurs with median impact pa- 
rameters for p°, u, (j) and J/ip production of 46 fm, 46 fm, 38 fm and 23 fm respectively, 
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corresponding to < 4 MeV/c for the light mesons and p± < 10 MeV/c for the J/ip. Below 
these pj_ values, interference is large. With lighter nuclei, the average impact parameters are 
slightly smaller, leading to higher momentum cutoffs. At the LHC, impact parameters are 
much larger, 170-290 fm for light vector mesons and 44-68 fm for the J/ip, and the reduction 
is large only for quite small values of p±. 

Despite the dramatically different p± spectrum, the the overall rate at ?/ = is unaffected 
by the interference. This is because {p±_) > {h/b), so the cos[p ■ b] oscillations average out, 
leaving the rate almost unaffected. Since the p± spectrum is almost independent of y, this 
remains true at other rapidities, so the total cross section is also unchanged. 

Figure 3 shows d?N/dp\ at ?/ = 1 for the same three cases. The magnitude of the 
interference is reduced because of the different amplitudes for the two directions. Besides 
the overall amplitudes, the spectrum is slightly different, because the photons have different 
energies, and so the p_[_ spectra are slightly different. 

This interference can be studied with the decays p° tt+tt^, uj tt+tt^, cf) K^K~ 
and J/ip ^ e^e^ which are readily reconstructible. Because of the high rates, the back- 
grounds to these processes should be small; simulations indicate that this is indeed the case, 
at least for a large acceptance detector |]T2[. 

There are a few details that may slightly reduce the interference. Because the photon 
emission is electromagnetic, it is possible for the emitting nucleus to be excited into a 
giant dipole resonance (GDR); since the Pomeron is uncharged, excitation is unlikely for the 
scattering nucleus. However, even for the photon emitting nucleus, the excitation probability 
is relatively low. Excitation is much more likely if the nuclei exchange an additional photon 
between themselves, with equal excitation probabilities for both nuclei, and not affecting 



the interference |]T3[. Radiative corrections and other higher order processes could also 
reduce the interference. These factors should be small because the nuclear trajectories are 
barely affected by the interaction. Likewise, backgrounds from non-resonant (but coherent) 
photonuclear and two-photon processes should be small. 

However, these backgrounds may introduce additional interference terms which may 



complicate the picture. For example, the process 77 — > e^e~ can interfere with e^e" from 
J/i}) decays; the rate for the 77 e+e" to mimic J/V' decays is small, but perhaps visible. 

In conclusion, interference between emission between two colliding nuclei dramatically 
changes the perpendicular momentum spectrum of exclusively produced vector mesons, sup- 
pressing production of vector mesons at low p^. Because the vector mesons decay before 
their wave functions can overlap, the interference involves all of the decay products, which 
never overlap with each other. So, the decay products have an entangled, non-local wave 
function. 

We would like to acknowledge useful conversations with Stan Brodsky and J0rgen Ran- 
drup. This work was supported by the US DOE, under contract DE-AC-03-76SF00098. 
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FIGURES 
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FIG. 1. Perpendicular momentum spectra for photons (dotted curves), Pomerons (dashed 
curves) and the final state vector mesons (solid curves) at (a) y = 0, and (b) y = —2 (corre- 
sponding to A; = 69 MeV in the lab frame) for (f) production in gold collisions at RHIC. The curves 
are each normalized to a maximum dN/dp\ = 1. Clear diffraction minima appear in the Pomeron 
spectra. Since k±^ is small compared to p_\_, the minima remain visible in the (j) p± spectrum. 
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FIG. 2. Expected p_L spectra for reconstructed cp and J/tp mesons at y=0 with (a,d) gold beams 
at RHIC, (b,e) silicon beams at RHIC and (c,f) calcium beams at LHC. The solid histograms include 
interference, while the dotted ones do not. Because of the smaller impact parameters production 
is peaked at higher pj_ in (b,e) than in (a,d). Because of the smaller impact parameters in (b,e), 
the interference dip extends to higher p± than in (a,d). In (c,f), the energies are higher, leading to 
higher impact parameters, pushing the interference dip to lower p±. The figure is normalized so 
that, without interference dn/dp\ = 1 at p_L = 0; the rates are given in Ref. ||l|. 
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FIG. 3. Expected p_L spectra for reconstructed cp and J/tp mesons at y=l with (a,d) gold beams 
at RHIC, (b,e) silicon beams at RHIC and (c,f) calcium beams at LHC. The solid histograms include 
interference, while the dotted ones do not. Because of the difference in amplitudes, dN/dp\ 7^ at 
Px_ = 0. The size of the dip at = depends on the ratio of amplitudes at y = 1 and y = —1; it 
is larger for lighter nuclei and higher collision energies. This figure is normalized so that, without 
interference dn/dp\ = 1 at p_L = 0; the rates are given in Ref. [ffl]. 
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